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Abstract 



We discuss properties of random fractals by means of a set of numbers that charac- 
fT) . terize their universal properties. This set is the generalized singularity spectrum that 

consists of the usual spectrum of multifractal dimensions and the associated com- 
plex analogs. Furthermore, non-universal properties are recovered from the study 
of a series of functions which are generalizations of the so-called energy integral. 
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A number of the properties of fractals are associated with their Hausdorff, 
box-counting, or fractal dimensions. But further information of a universal 
character is also encoded in secondary (singular) dimensions, some of which 
may be complex. For example, a recurrent theme in the study of fractals is 
that of asymptotic or logarithmic periodicity [1,2]. Most fractal objects and 
mathematical constructions thereof are not exactly scale invariant. Rather, 
they obey simple recurrence relations that relate an infinite but discrete set of 
scales. Recent physical examples include the appearance of complex exponents 
in diffusion-limited aggregation [3], crack propagation in two dimensions [4], 
and Boolean delay equations in the modelling of climate dynamics [5]. 

For instance, instead of the full scale invariance of a function of local variables, 

F{x) = X D F(x/X), (1) 
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we have the logarithmic analog of Bloch's theorem, 



F(x) = X F(x/X) G(m A), 



(2) 



where G may be a periodic function. 

This is a rather well-known relation but we feel that it has been under ex- 
ploited. The period of G is independent of the scaling dimension, D, and gives 
further information of the properties of the fractal object. 

The ways in which the fractal dimensions of the object manifest itself are 
manifold. Consider, in particular, a set of analytical quantities calculated from 
the real space distribution of a fractal object, namely, the set of correlation 
integrals: 



where q > 2, 9 is the Heaviside step function, and djj, is the measure of the 
object in question. The scaling properties of these correlation integrals with 
respect to the distance r define a countable set of dimensions forming part of 
the multifractal dimension spectrum [6-8]. 

The study of the scaling properties of such correlations is facilitated by con- 
sidering the corresponding energy integrals: 



where r is restricted to be less than the spatial dimension, d. Note that I q (r) 
is related to the Mellin transform of C q (r) [9]. 

At this point the spatial information of the fractal has been encoded in these 
energy integrals, which are typically, but not always, meromorphic. For partic- 
ular cases of deterministic Cantor sets, it is shown in [10,11] that the complex 
structure of these functions (4) reveal singularities that correspond to the rel- 
evant scaling dimensions of the theory. One usually one keeps only the most 
relevant, i.e., the one with the smallest real part, but the rest of the spectrum 
is important in studying finite-size effects. 

More precisely, the most relevant singularity is a pole on the real axis and 
has a numerical value that is a lower bound to the Hausdorff dimension. This 
was first proved by Frostman for q = 2 [12] (see also Falconer [13]). In some 
cases, the Hausdorff dimension also corresponds to the box-counting dimension 






(4) 



[13,14]. 
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Typically, the rest of the singularities are also poles, and appear as pairs of 
complex conjugates with real parts not smaller than the Hausdorff dimension. 
The imaginary parts of these poles correspond obviously to the logarithmic 
wavelength of the fractal, while the residues appear as the amplitudes of oscil- 
lations observed in the asymptotic scaling of various correlation integrals (3). 
This program of the analysis of a fractal object has been carried out, albeit 
in a somewhat scattered way, for the middle-third Cantor set and some of 
its deterministic generalizations [10,11,15]. This complex singularity spectrum 
has been called the multilacunarity spectrum [11]. It is the goal of our work 
to show that it is well-defined for classes of random fractals, and we explicitly 
compute the lacunarity of a particular example. 

It turns out that for some simple but important examples the Hausdorff di- 
mension is easy to calculate with the use of a little ingenuity. Just as simply, 
the satellite dimensions can be calculated in the same way without resorting to 
the explicit computation of the correlation integrals (3) or the energy integrals 

(4) . 

For the middle-third Cantor set, and for many other objects with simple re- 
cursive descriptions, we consider an equation that relates the relative scales, 
£i, and the relative (normalized) measures, p iy at successive levels of approxi- 
mation, 

i C i 

In the case of the middle-third Cantor set, ^=1,2 = 1/3 and Pi=i,2 = 1/2, and 

(5) becomes 2 9 ~ 1 /3 r = 1. The unique real solution, r(g) — (q — I)ln2/ln3, 
gives the Hausdorff dimension. In general, r = r(q) is not linear in q — 1 
and gives one the desired multifractal dimension spectrum via the relation 
D q = r(q)/(q — 1). In this way, the universal properties of the fractal, namely, 
its generalized dimensions, are rather easily computed [7,16]. 

However, as noticed in [11], a study of the middle-third Cantor set (and some 
deterministic generalizations), (5) also has complex roots. For instance, in the 
case of the middle-third Cantor set, 

In 2 271-1 

r{q) = (9 " I)--, + ^ J = 0. ±1. ±2, . . • (6) 

The imaginary parts of these complex roots correspond to the period of G in 
(2) and the logarithmic period observed in the correlation integrals [2]. 

The reason for the surprising success of this approach, which reduces the some- 
times formidable calculation of the energy integrals (4) to the computation of 
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a partition function (5), is that one has implicitly utilized the recursive struc- 
ture of the fractal distribution [encoded as (2)]. In doing so, one successfully 
captures the fact that the fractal has a very well defined set of real-space sin- 
gularities, and equates the determination of the spectrum of the many-body 
problem to the solution of a relatively simple transcendental equation. 

So far we have examined results for deterministic fractals, where the scale 
invariance of (2) is exactly satisfied. Consider now the case in which the fractal 
is not exactly self-similar, but is only statistically self-similar, i.e., scaling 
functions obey (2) only on average. We shall make precise what this averaging 
procedure entails (for related issues of averaging of stochastic hierarchical 
processes, see [17]). 

For a process that generates a generalized Cantor set by replacing each seg- 
ment at level / by to segments at level / + 1, we let the length of the seg- 
ments, £i(i = 1, . . . , to), be random variables with random probability mea- 
sures pi(i = 1, . . . , to). For this and other random fractals, we define a new 
set of correlation integrals as the expectations of (3), given a probability dis- 
tribution of the pj's and the ^'s: 

C q {r) = E \C q {r)} , (7) 

where E denotes expectation and C q (r) is the value of the correlation integral 
for a single realization of the random fractal. A new energy integral may be 
defined in precisely the same fashion: 

l q (r) = E[I q (r)}, r<d. (8) 

It has been shown by Falconer [18] that the probabilistic version of equation (5) 
still gives the relevant dimension spectrum, i.e, one has to solve the expectation 
equation, 

= 1, (9) 



to obtain the multifractal dimension spectrum. 

The conditions for the existence of a unique and meaningful real solution of 
(9) have been studied [18-20], and simple extensions of these considerations 
lead to the existence of well-defined complex solutions in complete analogy 
with the deterministic case. 

Consider then the following example of a randomized Cantor set. At level I of 
the recursive construction, we divide each segment into n equal segments and 



E 



E 



g 
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pick m of them at random with uniform probability. We assign to each of the 
m smaller segments a measure 1/m-th of the original segment. To simplify 
the presentation, we examine only the case for q = 2 and do not consider 
the more general model of Falconer [18], which involves possibly non-uniform 
probability distributions pi and £j. However, the computation may be gener- 
alized for higher-order correlation integrals and for non-uniform probability 
distributions satisfying the restrictions outlined in [18]. 

At level I of this process, the energy integral is simply related to the energy 
integral of the previous level: 

j(0( r ) =mn- T l (l - 1) (T) +TZ(m,n,T), (10) 



where the superscripts denote the finite-level approximation of the energy 
integral. The function 1Z is given by 

m , x 2n(m- l)(l-n T - 2 ) 

n(m, n, r = — ^ — — J —. 11 

V ' ' ; mn-1 r-1 r-2 v ; 



Explicit derivation of 1Z for this example is given as an appendix. Note that 
in (10), the prefactor mrC T is the expected value of the partition function for 
this model. So that the energy integral of the limiting distribution is simply 

J(r) = nm ^ T) . (12) 
w 1 - mn- T v 1 



All of the singularities of X(r) are given by the zeros of the denominator on 
the right-hand side of (12), 

mn~ T = 1, (13) 



as expected. It is easily checked that (11) is not singular at r = 1. 
The roots of (13) are at 

lnm .2irj . n , „ , ^ 

Tj = - + j = 0,±l,±2,... 

In n In n 



(14) 



Thus, we expect the correlation integral to exhibit oscillations of period In n. 
By using properties of inverse Mellin transforms [9], the correlation integral 
can be written as 



C(r)=r J, p + 2X; 



To 



cos 2nj 



lnr 
Inn 



(15) 
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Fig. 1. Log-log plot of r~ D C(r) versus r. Dotted curve is an ensemble average; solid 
curve is the expectation. 

where D = In mj Inn is the (second-order) correlation dimension, aj and <pj 
are real and are determined by the residues at Tj (for j non-negative): The 
residue at the j-th pole is of the form, aj exp(i(pj). We propose to call the 
modulus of aje^/rj (and the higher-order analogs) the lacunary amplitudes. 

In Fig. 1, we compare the expected scaling of the correlation integral with 
an ensemble average of the scaling for the case n = 3 and m — 2. The av- 
erage is performed over several (in this case, twelve) numerical realizations 
of the random fractal (approximated at level / = 15). We plot the residuals, 
r~ lnm / lnra C(r), versus r. The dashed line is the ensemble average, and the solid 
line exhibits the expected oscillations of (15). 

While it may seem superfluous to obtain the correlation integral from numer- 
ical computations once it has been calculated analytically, this exercise was 
interesting since we have not performed the average over a large ensemble. 
Rather we took the spatial average of a few instances of a random process. 
That both results were essentially the same over a number of logarithmic pe- 
riods is a rather natural self-averaging property of many fractal objects. The 
fluctuation of the ensemble average about the expectation may be described 
by higher-order correlation integrals and will be discussed in a forthcoming 
publication [21]. In the present case, the small size of the ensemble produces 
disagreement with predictions at r of the order of the system size. Further- 
more, numerical resolution affects the correlation for r smaller than e~ 13 . 
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We note that the measurement of the lacunarity spectrum is accessible using 
a variety of correlations. For instance, we may identify these complex dimen- 
sions in the logarithmically periodic oscillation of the Fourier transform [22] 
and the diffraction spectrum [22,23]. However, while the same logarithmic pe- 
riods are observed, the shape and phase of the oscillations vary. The periods 
arise from the additional discrete scale invariance of the underlying model and 
the numerical value of these periods are determined by the solutions of (9). 
The different measures of correlation (Fourier or Mellin transforms) reveal 
differently- valued residues located at the roots of (9). 

Previous studies of the inverse fractal problem, i.e., the extraction of the (pos- 
sibly stochastic) hierarchical process from the observed multifractal dimension 
spectrum, revealed ambiguities in the standard procedure: Namely, many mod- 
els can be made to fit a given multifractal dimension spectrum [16,24]. The 
approach described above provides the maximal characterization of the un- 
derlying multiplicative process without additional dynamical information. In 
this way, we may further distinguish between different hierarchical processes 
that give rise to fractals with similar dimension spectra (see also [17]). 

We have to stress that the lacunarity spectrum does not resolve all the inher- 
ent ambiguities. Many models, random or deterministic, can be made to fit a 
given lacunarity spectrum. However, using the lacunarity spectrum, we may 
distinguish between processes that have the same multifractal dimension spec- 
trum. What the lacunarity spectrum reveals is the possible additional discrete 
scale invariance which are not furnished by the previous attempts to char- 
acterize fractal systems. Furthermore, non-universal information is recovered 
from studying correlations (3) and energies (4). 

How does this work bear on fractal sets generated by low-order (deterministic 
or stochastic) dynamical systems? Preliminary work [21] suggests that the 
complex solutions of the Lyapunov partition function (a dynamical analog of 
(5) and (9), see [25]) describes the anomalous scaling observed in simulations 
(see, for instance, [1]). In addition, the fact that the lacunarity spectrum can 
be calculated from a partition function (5) immediately implies that periodic 
orbit expansions [26] can be used to calculate the spectrum. 

As for fractals generated by systems governed by large numbers of degrees 
of freedom (as featured most prominently in phenomena modeled by diffu- 
sion limited aggregation and in inhomogeneities of highly turbulent flows), 
this generalized multifractal description complements the traditional views. 
Studies thus far have taken the position that the deviations from strict power- 
law scaling are anomalous and, hence, have focused on establishing possible 
causes for this apparent deviation (an example being inertial range intermit- 
tency [8,27] in turbulent fluids). In our approach, log-periodic deviations may 
be accommodated rather naturally (see also [2,28]). Future efforts will be di- 
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rected towards the description of turbulent intermittency using the analysis 
followed in this paper. 

We gratefully acknowledge J. Fournier, R. Rosner, A. Sornborger, and E. 
Spiegel for helpful conversations. We thank R. Ball for bringing his work with 
R. Blumenfeld to our attention. We also wish to thank the anonymous referee 
for useful suggestions. This work was completed while F. J. S. was a Rosen- 
baum fellow at the Newton Institute. L. T. is supported by the U. K. Particle 
Physics and Astronomy Research Council. 

Appendix 

In this Appendix we present the details leading to the explicit formulas of the 
energy integral as given by (11) and (12). 

We need to consider first the behavior of the measure upon averaging. Since 
we have assigned an equal probability to every possible case of segmentation, 
the expectation value for the density, p(x)dx = dp(x), is uniform, i.e., 

E[p(x)] = l. (Al) 



To be able to perform the required multiple integrals, we need to evaluate 
the expectation of products of densities. Sufficient information about the joint 
distribution of these densities is obtained by considering one step in the re- 
cursive construction. After one such step, the interval L° = [0, 1] is divided 
into n subintervals Lj, % = 1, 2, . . . , n from which m subintervals will be chosen 
randomly. The set of all points (x, y) that belong to the subintervals Lj and 
Lj, respectively, form a square of size n _1 x n _1 in the x-y plane. Such a square 
will be labelled as shown in Figure 2. The shaded region R corresponds 

to those cases in which x and y are in disjoint intervals, i.e., i ^ j. 

Consider now one of the shaded squares, say The probability that the 

intervals Lj and L] are indeed chosen is simply j) / (m) ■ m ^ s case > eacn 
of the intervals will support a measure of total mass 1/m. Furthermore, since 
segmentation for each of these intervals proceeds uncorrelated, we have 



fn-2\ 

E [p(x)p(y)\ = ^g-Ei [p(x)] E 3 [p(y)} 



G) 

m{m — 1) n 2 
n(n — 1) m 2 ' 



(x,y)e(i,j), (A2) 



where E k denote expectation conditioned to the event that segment k is indeed 
chosen. 
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Fig. 2. Different regions of integration for the energy integral (10). Shaded regions 
(labelled by R) contribute to the function TZ(m,n,T); white regions contribute to 
X^(t). Labels of individual squares, are explained in the text. Shown here is 

a case for n = 5. 

Next, we consider the diagonal squares. These squares will contribute to the 
energy integral with probability (^-i)/(m) = m / n - The process of segmen- 
tation for a subinterval is identical to that of the original interval. Therefore, 
if L\ is chosen, the average of the product of densities p(x)p(y) restricted to 
this interval is identical to that of the original interval up to rescaling. The 
rescaling matches x from the L\ segment with the point x' = nx — i — 1 in L°. 
We have, for x, y G Lj, 

E [p(x)p(y)\ = (m/n)E t [p(x)p(y)\ 

= (m/n)(l/mfEW)p(y')} (A3) 

Note also that \x — y\ = (l/n)\x' — y'\. Thus the contribution to the energy 
integral from each of the diagonal squares is proportional to the overall ex- 
pectation of the energy integral 

= (m/n)(l/m) 2 n T l 2 (r) (A4) 



Summing over all squares we obtain relation (10) where TZ(m,n,T) can now 
be identified with the expectation value of the energy integral restricted to 
the shaded region. 



E 



dp(x)dp(y) 



\x 



[■<■:<) 
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To simplify the calculation of lZ(m, n, r) we note that 



E 



LR. 



dfi(x)dfi(y) 
\x — y\ T 



m(m 



n(n — I) 



R 



dxdy 
\x — y\ T 



(A5) 



Furthermore, the last integral satifies 

dxdy r dxdy x r f dxdy 



r r dxdy = r dxdy ^ r r dxdy 
J J \x — y\ T J \x — y\ T J J \x — y\ T 

1 R 1 1 00 1 1 



which readily gives the final result 



2nm-l ){l-n T - 2 ) 

R(m, n, r) = — - — . . n — J —. (A7) 

v ; m (n- l)(r- l)(r-2) v 7 
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